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Abstract 

The analytic asymptotic expressions for the Casimir free energy, pressure and en- 
tropy at low temperature in the configuration of one metal and one dielectric plate 
are obtained. For this purpose we develop the perturbation theory in a small pa- 
rameter proportional to the product of the separation between the plates and the 
temperature. This is done using both the simplified model of an ideal metal and of a 
dielectric with constant dielectric permittivity and for the realistic case of the metal 
and dielectric with frequency-dependent dielectric permittivities. The analytic ex- 
pressions for all related physical quantities at high temperature are also provided. 
The obtained analytic results are compared with numerical computations and good 
agreement is found. We demonstrate for the first time that the Lifshitz theory, when 
applied to the configuration of metal-dielectric, satisfies the requirements of ther- 
modynamics if the static dielectric permittivity of a dielectric plate is finite. If it is 
infinitely large, the Lifshitz formula is shown to violate the Nernst heat theorem. 
The implications of these results for the thermal quantum field theory in Matsubara 
formulation and for the recent measurements of the Casimir force between metal 
and semiconductor surfaces are discussed. 
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1. Introduction 



The Casimir effect [1] is the direct manifectation of zero-point oscillations of 
quantized fields. It finds multidisciplinary applications in quantum field the- 
ory, gravitation and cosmology, atomic physics, condensed matter and, most 
recently, in nanotechnology (see, e.g., the monographs [2,3,4,5] and reviews 
[6,7,8,9]). According to Casimir's prediction, the existence of zero-point oscil- 
lations leads to the polarization of vacuum in quantization volumes restricted 
by material boundaries and in spaces with non-Euclidean topology. This is 
accompanied by forces acting on the boundary surfaces (the so called Casimir 
force). The Casimir force acts between electrically neutral closely spaced sur- 
faces. It is a pure quantum phenomenon (there is no such a force in the frame- 
work of classical electrodynamics) being the generalization of the well known 
van der Waals force for the case of relatively large separations where relativis- 
tic effects become essential. 

The theoretical basis for the description of both the van der Waals and 
Casimir forces is given by the Lifshitz theory [10,11,12]. The main formu- 
las of the Lifshitz theory express the free energy and pressure of the van 
der Waals and Casimir interaction between two plane parallel plates as some 
functionals of the frequency-dependent dielectric permittivities of plate ma- 
terials. These formulas can be derived in many different theoretical schemes 
[2,3,4,5,6,7,8,9,10,11,12]. In particular, they were obtained in the framework 
of thermal quantum field theory in the Matsubara formulation [8]. During 
the last few years the Lifshitz theory was successfully applied to the interpre- 
tation of many measurements of the Casimir force between metal surfaces 
[13,14,15,16,17,18,19,20,21,22,23,24] and between metal and semiconductor 
[25,26,27,28]. 

A complicated problem of the Lifshitz theory is how to describe the Casimir 
interaction between real metals at nonzero temperature. The most convenient 
form of the Lifshitz formulas exploits the dielectric permittivity along the 
imaginary frequency axis. The latter is obtained from the tabulated optical 
data for the complex index of refraction by means of the Kramers-Kronig re- 
lations. The available data are, however, insufficient and must be extrapolated 
in some way to lower frequencies. In this respect the contribution from the 
zero frequency is of most concern. The point is that in Matsubara thermal field 
theory the zero-frequency term becomes dominant at large separations (high 
temperatures) whereas the contributions from all other Matsubara frequencies 
being exponentially small. In [29,30,31] the zero-frequency term of the Lifshitz 
formula was obtained by using the dielectric function of the Drude model. This 
results in a violation of the Nernst heat theorem in the case of perfect crystal 
lattices [32,33] and is in contradiction with experiments at separations below 
Ijim. [22,23,24,34]. The asymptotic value of the Casimir force at large sepa- 



2 



rations predicted in [29,30,31] is equal to one half of the value predicted for 
ideal metals, i.e., to one half of the so called classical limit [35,36]. 

Another approach [37,38] uses the dielectric permittivity of the plasma model 
to determine the zero-frequency term of the Lifshitz formula. This approach 
was shown to be in agreement with thermodynamics [32,33] and consistent 
with experiment [23]. It predicts the magnitudes of the Casimir force at short 
separations in qualitative agreement with the case of ideal metals. At large 
separations the predicted force magnitude is practically equal to that for ideal 
metals. Very similar results, which are also in agreement with the require- 
ments of thermodynamics and consistent with experiment, are predicted by 
the surface impedance approach [39,40]. The controversies among different ap- 
proaches to the thermal Casimir force between metals are detailly discussed 
in [31,33,40,41,42,43,44]. 

Recently it was demonstrated [45,46,47] that even the traditional application 
of the Lifshitz formula to the case of two dielectric semispaces presents prob- 
lems. In [45,46,47] the analytic asymptotic expressions for the Casimir free en- 
ergy, pressure and entropy at low temperatures (short separations) were found 
for two dielectrics. It was shown that if the dielectric materials possess finite 
static dielectric permittivities the theory is self-consistent and in agreement 
with thermodynamics. If, however, a nonzero dc conductivity of dielectrics 
is taken into account (any dielectric at nonzero temperature is characterized 
by some nonzero dc conductivity which is many orders of magnitude lower 
than for metals) , this leads to a qualitative enhancement of the Casimir force 
and a simultaneous violation of the Nernst heat theorem. (Note that the dc 
conductivity of dielectrics was taken into account in [48] to explain the large 
observed effect in noncontact friction [49].) In [45,46,47] the phenomenological 
prescription was proposed that the dc conductivity of dielectrics is not related 
to the Casimir interaction, and to avoid contradictions with thermodynamics 
it should not be included in the model of dielectric response. 

The difficulties which were met in the application of the Lifshitz theory to two 
metal and two dielectric plates attracted attention to the case when one plate 
is metallic and another one dielectric. This configuration was first considered 
in [50]. It presents the interesting opportunity to investigate the Casimir force 
in the case when different plates are described by quite different models of the 
dielectric response. In [50], however, only the first leading terms in the low- 
temperature asymptotic expressions for the Casimir free energy and entropy 
were obtained, and the pressure was derived only in the dilute approximation. 
In the analytical derivations in [50] (see also the review [47]) it was supposed 
that the metallic plate is made of ideal metal and the dielectric of the other 
plate is described by the frequency independent dielectric permittivity. These 
suppositions narrow the applicability of the obtained results. Also, the role of 
the dc conductivity of a dielectric plate was not investigated for plates with 
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frequency-dependent dielectric permittivities. 

In the present paper we develop the analytic approach to the thermal Casimir 
force acting between metal and dielectric permitting to find several expansion 
terms in the asymptotic expressions for all physical quantities at low temper- 
ature. This approach is applied not only to the configuration of ideal metal 
and dielectric with frequency independent dielectric permittivity but also to 
real metal and dielectric described by the dielectric permittivities depending 
on frequency. We pioneer in derivation of the low-temperature asymptotic ex- 
pressions for the Casimir free energy, entropy and pressure between real metal 
and dielectric. The asymptotic behavior of all physical quantities at high tem- 
peratures (large separations) is also provided. What is more, the obtained rep- 
resentation for the Casimir free energy permits to find the low-temperature 
behavior of the Casimir force acting between a metal sphere and a dielectric 
plate (or, alternatively, dielectric sphere above a metal plate). This can be 
done with the help of the proximity force theorem [51]. The configuration of a 
sphere above a plate is most topical in experiments on the measurement of the 
Casimir force [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. (Note that for 
the experimental parameters the error introduced by the use of the proximity 
force theorem was recently shown to be less than 0.1% [52,53,54,55].) Thus, 
our results will find immediate utility in experiment. The analytic expressions 
for the Casimir interaction between metal and dielectric at zero temperature 
are also found here for the first time. We determine the applicability region of 
the obtained analytic formulas by compairing them with numerical computa- 
tions using the tabulated optical data for metaUic and dielectric materials. The 
fundamental conclusion following from our results is that the Lifshitz theory, 
applied to the configuration of a metal and a dielectric plate, is in agreement 
with the Nernst heat theorem if the static dielectric permittivity of a dielectric 
plate is finite. Note that this conclusion could not be achieved by using the 
numerical computations which inevitably identify zero with all nonzero num- 
bers in the limits of a computational error. If, however, the dc conductivity of 
a dielectric plate is included in the model of dielectric response, we show that 
the Nernst heat theorem is violated. This is in analogy to the same conclusion 
in [45] obtained for the configuration of two dielectric plates and confirms our 
phenomenological prescription that the dc conductivity is not related to the 
van der Waals and Casimir forces and should not be included in the model of 
dielectric response. Recently this prescription was confirmed experimentally 
[28]. 

The paper is organized as follows. In Section 2 we summarize the main formu- 
las of the Lifshitz theory for the configuration of one plate made of metal and 
another one made of dielectric. Section 3 is devoted to the simplified model 
where the metal is an ideal one and dielectric is described by a constant di- 
electric permittivity. In the framework of this model a perturbation formalism 
applicable at low temperatures (short separations) is developed. In Section 4 
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the realistic case is considered when the dielectric permittivities of both metal 
and dielectric plates depend on the frequency. The analytic asymptotic ex- 
pressions for the free energy, entropy and pressure of the Casimir interaction 
at both low and high temperatures are obtained. Section 5 contains the com- 
parison between the analytical results and numerical computations using the 
tabulated optical data for plate materials. The application region of the de- 
rived asymptotic expressions is determined. In Section 6 it is shown that the 
inclusion of the dc conductivity in the description of dielectric plate leads to a 
violation of the Nernst heat theorem. Section 7 contains our conclusions and 
discussion. 



2. Lifshitz formula in the configuration of metal and dielectric 
plates 



We consider two thick parallel plates (semispaces) at temperature T in ther- 
mal equilibrium separated by the empty gap of width a. One plate is made 
of metal with the dielectric permittivity e^{uj) and another of dielectric with 
permittivity e^iuj). The free energy of the van der Waals and Casimir in- 
teraction between the plates per unit area is given by the Lifshitz formula 
[10,11,12,45,50] 



I. rp OO 



271 



1=0 



k\dk\ 



x{\n[l-rf{iiM)r^{iiM)e-^'^'] 
+ In [l - rf (6, Mr^^te, Me-'^ } • 



Here the plates are perpendicular to the z axis, k^^ — \k±\ is the magnitude of 
the wave vector in the plane of plates, = 2nkBTl/h are the Matsubara fre- 
quencies, and ks is the Boltzmann constant, '''^i^ are the reflection coefficients 
for metal (M) and dielectric (D) plates for the two independent polarizations 
of electromagnetic field calculated along the imaginary frequency axis. Index 
II stands for the electric field parallel to the plane formed by k± and the z 
axis (transverse magnetic field), and index _L stands for the electric field per- 
pendicular to this plane (transverse electric field). The explicit expressions for 
the reflection coefficients are [45,50] 



M,D 7 M.D 

£i ' qi-ki ' 
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where 



Qi 



M,D 



(3) 



and 



e, — e 



M,Di 



(4) 



The pressure of the van der Waals and Casimir interaction between metal and 
dielectric (i.e., the force per unit area of plates) is obtained from 



Using Eq. (1) we arrive at 



P(a,r) = fl - ^^0^) lki.dk^qi (6) 



Using the proximity force theorem [51], one can obtain from Eq. (1) the ap- 
proximate expression for the Casimir force acting between a sphere and a 
plate 

F(a,T) = 27rRJ^{a,T). (7) 



This equation is widely used for the interpretation of measurements of the 
Casimir force [13,14,15,16,17,18,19,20,21,22,23,24,25,26]. Recently both exact 
analytic and numerical results for the Casimir force in the configuration of 
a cylinder above a plate (electromagnetic case) and for a sphere above a 
plate (scalar case) were obtained [52,53,54,55]. It was shown that the error 
introduced by the use of Eq. (7) for the experimental parameters in already 
performed experiments is less than 0.1%. Using Eq. (7), the analytical results 
derived below for metal and dielectric plates, can be immediately applied for 
the interpretation of measurements of the Casimir force between Au coated 
sphere and Si plate [25,26,27,28]. 
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The analytic perturbation expansions in Eqs. (1) and (6) can be conveniently 
performed by using the dimensionless variables Q and y 



= — = — = v = ^m, 

UJc C 



(8) 



where uJc = c/ (2a) is the characteristic frequency of the Casimir effect and 
T = A-KkBoT /{tic). In terms of these variables the free energy (1) takes the 
form 



+ ln 



l-rf(0,?/K(0,Z/)e-^]} 



Using the variables (8), the reflection coefficients (2) are 



(9) 



^11 iChy) - 



' y + \/y^ + Wi - 1) 



where in accordance with Eq. (4) sf^'^ = e^'^{iC,iUc). 
The pressure (6) is rearranged as follows: 



(10) 



P(a,T) 



y'^dy 



X 



rf(0,l/)rf(0,l/) , r^{Ci,y)r^{Ci,y) 



+ 



ey-rf{0,y)rf{Ci,y) - (0, |/)rf (0, |/) 



The other important characteristic of the van der Waals and Casimir interac- 
tion is the entropy 



S{a,T) 



or {a, T) 
dT 



(12) 



In [32,33] the behavior of the Casimir entropy at T — > was used as a phe- 
nomenological constraint on the selection of theoretically consistent models of 
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the dielectric response for real metals at low frequencies. It was proposed that 
all consistent models should satisfy the thermodynamic condition S{a, 0) = 0, 
i.e., be in agreement with the Nernst heat theorem. In [45] it was demonstrated 
that this condition is respected for two dielectric plates with the finite static 
dielectric permittivities. The new analytic expressions for the free energy ob- 
tained in the present paper permit investigate the behavior of entropy in the 
configuration of one metal and one dielectric plate and find when it vanishes 
with vanishing temperature. 



3. Model of ideal metal and dielectric with constant dielectric per- 
mittivity 



To find the analytic expressions for the free energy, pressure and entropy of 
the Casimir interaction between metal and dielectric, we start from a simpli- 
fied model when the metal is an ideal one and the dielectric possesses some 
finite dielectric permittivity Eq independent on the frequency. Such modeling 
is widely used in Casimir physics (see, e.g., [2,3,5,6,8,10,11,12,56]). It provides 
rather good description of real metals and dielectrics at sufficiently large sepa- 
rations between the interacting surfaces. For an ideal metal it holds \s^\ — oo 
at all frequencies and from Eq. (10) one obtains 

rr(0,y) = l, rf(0,y) = l, l>0. (13) 



Using Eq. (13), the free energy (9) and pressure (11) are represented in a more 
simple form. 



X {in [l - rf (0, y)e-y] + In [l - (0, y)e-'] } , 



X 



(14) 



Notice that in the framework of our model in Eq. (10) it holds = Eq, 
i.e., the dielectric permittivities computed at different imaginary Matzubara 
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frequencies do not depend on L In particular, at Z = it follows: 



(15) 



In fact Eq. (15) is valid not only for our simplified model but for any dielec- 
tric with a finite static dielectric permittivity £^(0) = Sq < oo. Usually for 
nonpolar dielectrics e^iiE) = Eq = const in the frequency region from ^ = 
up to rather high frequencies of about lO^^rad/s and for higher frequencies 
s^{i^) decreases to unity. The simplified model does not take the latter into 
account (in the next section we show that this does not influence the first 
terms in the asymptotic behavior of the free energy, entropy and pressure at 
low temperature). 

Now we proceed with the derivation of the asymptotic behavior of the Casimir 
free energy at low temperature (r <^ 1). Using the Abel-Plana formula [3,8] 



'=0 



- 1 



(16) 



where F{z) is an analytic function in the right-plane, we can rearrange Eq. (14) 
to the form 



:r(a,r) = £;(a) + Ajr(a,r). 



(17) 



Here, 



E{a) 



he 



oo oo 



327r2a3 



dC f{Cy)dy 



(18) 



is the energy of the Casimir interaction at zero temperature. 



Ajr(a, T) 



ihcT 



dt 



F{iTt) - F{-iTt) 



(19) 



is the thermal correction to it, and the following notations are introduced. 



/(C, y) = y In [l - rf (C, y)e-^] + y In [l - rf (C, y)e-^ 

oo 

F{x) = j dyf{x,y). 



(20) 
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The expansion of f{x, y) in powers of x takes the form 



f{x,y) = yln(l - roe ^) 



(21) 



X 



I - 1 



-I? oo 



4y 



^0^ + ln=l 



-ny 



y 



+ 0{x^). 



To find F{x) in Eq. (20) we integrate the right-hand side of Eq. (21) with 
respect to y. Notice that the first term on the right-hand side of Eq. (21) does 
not contribute to the first expansion orders of F{ix) — F{—ix) which is in fact 
the quantity of our interest. This is because in the expression 



j ydy]n{l - roe"^) = j vdvln{l - roe-"") + 0{x'^), 



(22) 



where the new variable v = y — x was introduced, the first-order in x con- 
tribution vanishes. Thus, this term could contribute to F{ix) — F{—ix) only 
starting from the third expansion order. Integrating the second term on the 
right-hand side of Eq. (21) using the formula 



dy = —Ei{—nx), 



(23) 



where Ei{z) is the exponential integral function, we finally obtain 

F{ix) - F{-ix) = m ^^^'^^ x" - i-fx^ + O(x^), 

4(£o + 1) 



(24) 



where the third order real expansion coefficient 7 cannot be determined at 
this stage of our calculations because all powers in the expansion of f{x, y) in 
powers of x contribute to it. 



Now we substitude Eq. (24) in Eq. (19) and find the free energy (17) 



T{a, T) = E{a) - 



tie 



327r2a3 



C(3) (gp 1)^ ^3 li' ^4 I o/'-^5\ 



(25) 



where — 7/240 and C,[z) is the Riemann zeta function. 

The Casimir pressure is obtained from Eqs. (5) and (25). It is equal to 



P(a,r)=Po(a)- 



he 



327r2a4 



X4t' + 0(t5) 



(26) 
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where Po{a) — —dE{a)/da. 



In order to determine the coefficients K4, we now start from the Lifsliitz 
representation of the pressure in Eq. (14). Using the Abel-Plana formula (16), 
we rearrange Eq. (14) to the form analogical to (17)-(19), 



P{a,T)^Po{a) + AP{a,T), 



Po{a) 



3hc 
32^ 



00 00 



d( f{C,y)dy, 



(27) 



AP{a,T) 



c 

ihcT J ^^^{irt) - ^{-irt) 



327r2a4 



^2nt _ 1 



Here Po{a) is the Casimir pressure at zero temperature, AP(a, T) is the ther- 
mal correction to it and the following notation is introduced: 



$11. 1 (x. 



00 



y'^ri±{x,y) 
ey - ri±{x,y)' 



(28) 



To find the expansion of ^{ix) — $(— ix) in powers of x, we first deal with 
$^(x). By adding and subtracting the asymptotic behavior of the intergrand 
function at small x, 



— r±(x, y) 4 



(29) 



under the integral in Eq. (28) and introducing the new variable v = y/x, the 
function $x {x) can be identically rearranged and expanded in powers of x as 
follows: 



00 

$^(x) = ^{£q - l)x^e-'^ + x^ Jdv 



E 

n=l 



^\{s^-l)x'{l-x) 



(30) 



+ x / dv 



v'^rj_{v) - 1 



1 - r±{v) 4 



+ 0(a;^). 



Calculating the integral on the right-hand side of Eq. (30), we arrive at the 
result 



$^(x) = '-l^x' - I ( e^^e^ -l]x^ + 0{x^). 



(31) 
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To deal with ^\\{x), we add and subtract Tinder the integral in Eq. (28) the 
two first expansion terms of the integrated function in powers of x, 



$11 (x) = j y^dy 

X 

oo 

+ / V^dy 



r\\{x,y) ro 



+ 



e^roe-yx^ 



ey-r\\{x,y) ey-ro y^{e^ + - roe-yf 



.(32) 



The asymptotic expansions of the first and second integrals on the right-hand 
side of Eq. (32) are 



2Li3(ro) 



2(^0^ + 1) 



+ -{e^ - - 2)x' + Oix% 



- 1) - le^ie^^f^ - 1) + ^e^ie^ - 1)^ 

+ 0{x'), 



X 



(33) 



(34) 



respectively, where Lin{z) is the polylogarithm function. By summing Eqs. (33) 
and (34) we find 



$,Kx) = 2Li3(ro)- 2^^,^^^ 



(35) 



So" - 1) + - 1) - ^e^ie^ - l)V^o" 



x^ + 0{x^). 



After summing Eqs. (31) and (35) the following result is obtained: 

^ix) - ^-ix) = -| (l - '2e^^+ {e^y ^) + 0{x'). (36) 



Substituting this in Eq. (27) and integrating we arrive at the asymptotic 
expression for the Casimir pressure in the limit of small r, 



P{a,T) = Po{a) - 



he 



327r2a4 



360 



(37) 



The comparison of this equation with Eq. (26) leads to the explicit expression 
for the coefficient K^: 



(38) 
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and, thus, the exphcit asymptotic expression (25) for the free energy is also 
fully determined. 

Notice that the energy and pressure at zero temperature [E{a) and Po{a) 
defined in Eqs. (18) and (27), respectively] depend on the separation through 
only the factors and in front of the integrals. They can be conveniently 
presented in the form 



E{a) = -—-i^DM{e^). (39) 
^o(a) = -— ^V'dm(£o)> 
where the function iPdm{£o) is defined as 

oo oo 

^DM{e^) = / / /(C, y)dy. (40) 



In fact, ipDM in Eqs. (39), (40) is the correction factor to the famous Casimir 
result [1] obtained for two ideal metals. It is equal to the function (poM intro- 
duced in [12], multipled by vq. 

The function il^DMi^o) in Eq. (40) can be presented in a more simple analytical 
form as follows. Presenting the logarithms in Eq. (20) as series and changing 
the order of integrations, one obtains 



i/Jdm{£o ) 



45 
2^ 

y 



- / ydye- 

n=l n 



■ny 



(41) 



dC{[r^{C,y)] +[rf(C,?/)] }■ 





Introducing the new variable w — (^/y, we rearrange Eq. (41) to the form 



45 °° 1 7 

n=l n 



-ny 



(42) 



X 



1 

ldw{[r^{w)Y+[r^{w)Y}, 



where 
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-D 



1)^2 



(43) 



Calculating the integral in y and performing the summation with respect to 
n in Eq. (42) one arrives at 



i^DMiso) ^ J dw^Lu rf{w) +U4 r^{w) }. 



(44) 



In Fig. 1 the function (44) is plotted versus £^ as a solid line (when — > 1 it 
goes to zero and when £^ — > 00 it goes to unity reproducing the limit of ideal 
metals). 

It is notable that the model under consideration represents correctly the 
Casimir energy and pressure (39) at T = in only the retarded regime (i.e., at 
sufficiently large separations). As to the thermal corrections in Eqs. (25) and 
(37), the obtained expressions are valid also at short separations in a nonre- 
tarded regime under the condition that the parameter r is sufficiently small 
due to sufficiently low temperature. 

From Eqs. (12) and (25) the asymptotic behavior of the Casimir entropy in 
the limit of small r is given by 



S{a,T) 



1287r3a2 (e 



+ 1) 



(45) 



X 



1 - 



135C(3) (£? - if 



T + 0(t' 



As is seen from Eq. (45), the entropy of the Casimir interaction between 
metal and dielectric plates vanishes with vanishing temperature as is required 
by the Nernst heat theorem (note that the first term of order in Eq. (45) 
was obtained in [50]). The important property of the perturbation expansions 
in powers of r in Eqs. (25), (37) and (45) is that it is impermissible to consider 
the limiting case 5:^ ^ 00 in order to obtain the case of two ideal metals like 
it was discussed above in application to Eq. (39). The mathematical reason 
is that in the power expansion of functions depending on Eq as a parameter 
the limiting transitions Eq ^ 00 and r — > are not interchangeable. Of great 
importance is the possibility to apply Eq. (45) at as small T as is wished. This 
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is the principal advantage of analytical calculations as compared to numerical 
ones. 

Now we consider the opposite limiting case r ^ 1, i.e., the limit of high 
temperatures (large separations). Here the main contribution to the free en- 
ergy (14) is given by the term with I — whereas all terms with I > 1 are 
exponentially small [8], 





By integrating in Eq. (46) we obtain 



For the Casimir pressure and entropy at r ^ 1 from Eqs. (5), (12) and (47) 
it follows 



4. Thermal Casimir force between metal and dielectric with 
frequency-dependent dielectric permittivities 



In this section we obtain the analytic expressions for the low-temperature 
behavior of the Casimir interaction between metal and dielectric plates taking 
into account the dependence of their dielectric permittivities on the frequency. 
The metal plate is described by the dielectric permittivity of the plasma model, 



-M 



(49) 



where ujp — 27rc/Ap is the plasma frequency, and \p is the plasma wavelength. 
In the theory of the thermal Casimir force this description was first used in 

[37,38] and was shown to work good at separations between plates greater 
than the plasma wavelength. At such separations the characteristic frequency 
of the Casimir effect ujc belongs to the region of infrared optics where the 
relaxation processes do not play any role [57]. 
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For dielectric plate we use the Ninham-Parsegian representation of the dielec- 
tric permittivity along the imaginary frequency axis [58,59], 

e^m = ^+i:^- (50) 



Here Cj are the absorption strengths satisfying the condition 

= (51) 



and ujj are the characteristic absorption frequencies [recall that now Sq = 
£^{0) < oo]. Eq. (50) gives a very accurate approximate description of the 
dielectric properties for many dielectrics. It has been successfully used by 
many authors for the comparison of experimental data with theory [60] . 

From Eq. (50) we return to the same values (15) of the reflection coefficients 
of the dielectric plate at zero frequency as were obtained in the simplified 
model of the frequency-independent dielectric permittivity. Thus, due to the 
zero value of r^{0,y), the transverse electric mode at zero frequency does not 
contribute to the free energy (9) of the Casimir interaction between metal and 
dielectric regardless of the value of r^{0,y) for a metal. As was told in the 
Introduction, there are different approaches on how to correctly calculate the 
transverse electric coefficient at zero frequency, r^(0, y), for a plate made of 
real metal. In the configuration of metal and dielectric this problem, however, 
does not influence the result. Note that if we would use instead of Eq. (49) the 
Drude model, taking relaxation into account, the prime perturbation orders 
in all results below remain unchanged for metals with perfect crystal lattices. 
The role of impurities in the validity of the Nernst heat theorem in the case 
of two metal plates is discussed in [31,33,41,42,43,44,61]. 

We start from Eq. (9) for the free energy. Once again, using the Abel-Plana 
formula, Eq. (9) can be represented by Eq. (17) as the sum of E{a) in Eq. (18) 
and AjF(a, T) in Eqs. (19) and Eq. (20), where we mark by a hat all quantities 
related to real metal and dielectric. The single difference is that the function 
f{x,y) in Eq. (20) should be replaced by 



f{x,y) = f\\{x,y)+ f±{x,y), 



(52) 



I II I 



It is notable that for real metal and dielectric E{a) and A^{a,T) in Eq. (17) 
may lose the obvious meaning of the energy at zero temperature and the 
thermal correction to it. In fact, this meaning is preserved only in the case 
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when the dielectric permittivities e^'^(i^) do not depend on the temperature 
as a parameter hke it was in Section 3. In the latter case it holds 

Aj^(a, T) = T{a, T) - T{a, 0) = T{a, T) - E{a) (53) 



in accordance with the intuitive definition of the thermal correction. If, how- 
ever, e^^{i^) or s^{i^) or both depend on the temperature as a parameter, 
Eq. (53) is violated. In this case AjF(a,T) defined in Eq. (19) takes into ac- 
count only the part of temperature dependence of the free energy originating 
from the Matsubara frequencies and is not equal to ^{a, T) — T{a, 0). More- 
over, in this case E{a) in Eqs. (17) and (18) is in fact temperature-dependent 
and it would be more correct to use the notation E = E{a, T). 

To obtain the analj^ic expressions of our interest we develop the perturbation 
theory in two small parameters r and r] = S/{2a), where S = \p/(27i) is the 
penetration depth of the electromagnetic oscillations into a metal. For the sake 
of simplicity we will consider dielectrics which can be described by Eq. (50) 
with only one oscillator, i.e., with j = 1. The high-resistivity Si is a typical 
example of such materials. The function F{x) in Eq. (20) can be conveniently 
presented in the form 

oo 

F(x) = F\\(x) + F^(x), F||,^(x) = I dyf\\,^(x,y). (54) 

X 



As a first step we perform the expansion with respect to the powers of small 
parameter rj. This results in: 



oo oo 

F||(a;) = j ydy In 1 — f|j^(x, y)e~^ +2x^r] j dy 



f^{x, y) 
qv — fP{x, y) 



oo 



e^fji (x, y) 



y[ey — fF{x, y)]' 



oo oo 

F±{x) = j ydy In [l - ff (x,y)e~^] +"2^1 v'^dy 



r±ix,y) 



ey — fj^{x, y) 



oo 

277' / y^dy 



eyfj_{x, y) 
[ef-rf(x,y)]^ 



+ 0{rf). 



(55) 



The dielectric reflection coefficients in Eq. (55) are obtained after the substi- 
tution of Eqs. (50) and (51) with j = 1 in Eq. (10): 
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where b = uJc/oji. 



Let us consecutively consider the contributions to F{x) from the terms of 
order rf, rj and rj'^ in Eq. (55). As to the terms of order 77° [the first and fourth 
hnes in Eq. (55)], the expansion in powers of small x (small r) performed using 
Eq. (56) leads to 



F^o{x) ^ F{x) -b'x^ < 



nrg Ei(—na;) 



(57) 



n=l 



-rl Y: ^^oEi [-{n + m + ^-^Ei(-x) + 0{x'). 



n=l 



Here F(x) was already calculated in Sec. Ill and results in Eq. (24). The 
additional contributions to the right-hand side of Eq. (57) lead to the term of 
order in F^o(irt) —Fjf){—iTt) and of order in the free energy. Thus, they 
can be omitted (like in Sec. Ill, we preserve only the terms of order and 

To find the contribution to F{x) of order rj [we use the notation F^(a;)], we 
expand in powers of a; the following quantities under the integrals in Eq. (55): 



X X 



ey — fP{x,y) ey — Tq y'^{eQ -\- l){ey — TqY 



=y 



26V -^^ ^ + 0(x^), (58) 

(5? + l)(e^-ro)2 ^ ^' ^ ^ 



y'r^{x,y) _ , {e^ - l)e-y ,{e^ -\fe-''y{2ey -\) 



2 



— (x, 7/) 4 16|/ 

4 

By integrating of the third terms on the right-hand side of equations (58) 
with respect to y from x to infinity, we find that they contribute to Fy^{iT^) — 
Fj^{—iTt) only in the order and, thus, to the free energy in the order r^. 
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Because of this they can be omitted. The integration of the first two terms on 
the right-hand side of equations (58) leads to 



- = irirH^ie^ - l){e^ + 2), (59) 

Prom Eq. (59) it follows 

F„{irt) - F,{-irt) = irjrH' _ (60) 



As to the terms of order rj"^ in Eq. (55), their lowest order contributions to 

ii(irt) — ||(— irt) and to _|_(irt) — -F^2 ±(— irt) are of order r"^ and r^, 
respectively. This leads to the respective contributions of order and to 
the free energy which we omit in our analysis. 

Using Eq. (19), the respective correction to the Casimir free energy takes the 
form 

fir 

^^-^"'^^ = ~ 307207r%3 ^^'(^o^ " ^^^^'o + H)- (61) 



Remarkably, r^r^ ~ and the correction (61) does not depend on the sep- 
aration. Thus, there is no correction to the Casimir pressure of order r/r' 
with g < 4 due to the finite conductivity of a metal plate. Recall that in the 
configuration of two ideal metal plates the main thermal correction at low 
temperature is of order r*^. If the nonidcality of a metal is taken into account, 
the correction of order arises [8]. From this it follows that the thermal 
correction in the configuration metal-dielectric is less sensitive to the finite 
conductivity of a metal than in configuration of two metals. 

Combining the contributions from the zeroth and first orders in t] in Eqs. (25) 
and (61), the free energy at low temperatures for the configuration or real 
metal and real dielectric is 



T{a,T)^E{a) 



he 



327r2a3 



C(3) 



.D 



If 



167r2 ej? + l 



(62) 



+ 960^"^"^^^^'^ + ^^^'^^ 



where K4 is defined in Eq. (38). It is notable that the low-temperature behavior 
of the free energy is not influenced by the absorption bands of the dielectric 
material and are determined by only the static dielectric permittivity. This is 
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in analogy to the case of two dielectric plates [45,46]. For the Casimir pressure 
between plates made of real metal and dielectric Eq. (37) is preserved with 
the replacement of Po{a) for Po(a) given below. 

Now we derive the analytic representation for the Casimir energy E[a) in the 
configuration with one plate made of real metal and another plate made of 
real dielectric. Expanding in powers of rj in Eq. (18), with / replaced by / 
from Eq. (52), we obtain 



E(a) 



he 



327r2a3 
+ ln (l-ff(C,?/)e-^) 

oo oo 

+ 2nldc[ejdy-^ 

oo oo 

-27]' J dc[c'J dy 



dC / ydy [in (l - ff{C,y)e~ 



(63) 



+ 



rf(C,?/) 
'^-rf(C,2/) ^ 

e'rfiCy) 
y{ey-f^{C,y)f 



oo 

/ ^^^^( 



r^iCy) 



+ / y^dy 



fnCy) 



e'rIiCy) 



{ey-rfiC,y)y 



> . 



Here, the reflection coefficients for dielectric with the frequency-dependent 
dielectric permittivity are defined in Eq. (56). For many dielectrics, admitting 
the presentation (50) with one oscillator, the characteristic frequency at typical 
separations is much less than the absorption frequency leading to 6 = uJc/oji <^ 
1. In fact the small parameter h is of order of another small parameter r]. The 
expansion of Eq. (56) in powers of h takes the form 



rn (x, y) = rii [x, y) — o 



+ 0(6^), 

(64) 



f±{x,y) = r1{x,y) - 6^ 



{e^ - l)x'y 



- l)x^ + y'(y + ^{e^-l)x^ + y^ 



where r^j_{x,y) are the refiection coefficients for dielectric with a frequency- 
independent dielectric permittivity e^. Our goal is to obtain the expansion of 
£■(0) up to the second powers in the small parameters 77 and b. 

To attain this goal, we note that Eq. (64) contains the zeroth and second 
powers in b. Thus, both of them should be substituted in the zeroth power 
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in T) in Eq. (63). Considering the terms of order rj and r)"^ in Eq. (63) we 
should restrict ourselves by only zeroth order in b, i.e., replace r\\^±{x,y) for 
r\\^±{x,y). The calculational scheme of all coefficients accompanying rj, rf and 
6^ is the same as was used in Section 3 for obtaining the analytic expression 
for the function ipoMi^o)- It consists in the expansion of the integrands in a 
power scries, changing the order of integrals and introducing the new variable 
w = C,/y. In the order rf in Eq. (63) we obtain the contribution ah^cady 
calculated in Eqs. (39), (44) and the contribution of order 6^. The latter takes 
the form 



(65) 



X < 



2 + {e{^ -2)w^] \rP{w) 





n-l 



W 



+ 



n-l 



} . 



After the integration with respect to y and summation we arrive at 



Eh2[a) = . „ ^ / aw 



w 



X 



2 + 



(66) 



+ U5 



rf(w) 



Following the same procedure for the terms of order 77 in Eq. (63) we obtain 
1 



Eri{a) 



ShcT] 

8^ 



dw |w^Li4 (w) 



Lid 



n 



(67) 



Quite analogically for the terms of order 77^ in Eq. (63) it follows 

1 



Er,2{a) 



ShcT] 

'2^ 



dw ^w'^L 
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rli'iw) 



Lid 



n 



(68) 



By combining Eqs. (39) and (66)-(68) one arrives at the Casimir energy in 
the configuration of metal-dielectric plates made of real materials, 



E(a) 



77^ he 
'720^ 



D 



(69) 
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where the positive coefficients Ci, C2 and B are defined as 



Er,{a) 
■ E{a) ' 

" E{a) 



E{a) ' 



(70) 



and E{a) is given in Eq. (39). 



In Fig. 1 the above coefficients arc plotted as functions of Sq by the long- 
dashed lines 1 and 2 (Ci and C2, respectively) and by the short-dashed line 
[B). From Eq. (69) and Fig. 1 it is easy to obtain the respective analytic 
expression for the Casimir pressure, 



Po(a) = 



77^ he 

'240a4 



iPdm{£o) 



.(71) 



Equations (69) and (71) give the possibility to simply find the Casimir energy 
and pressure between metal and dielectric with rather high precision (see the 
next section). 

Prom Eqs. (12) and (62) we obtain the asymptotic behavior of the Casimir 
entropy at small r for metal-dielectric plates made of real materials, 



S{a,T) 



3A;bC(3)(£ 



D 



USiT^a^ie^ + 1; 



-t'<1 



45C(3)(£?-1) 



r 



(72) 



X 



2e 



D L,D 



V '^o 



3(^^ - 1) 



-{5e^ + ll)r} 



+ 0(r^) 



> . 



As is seen from Eq. (72), S{a, T) goes to zero when the temperature vanishes as 
is required by the Nernst heat theorem. This completes the proof of the impor- 
tant statement that the Lifshitz theory in the configuration metal-dielectric 
is consistent with thermodynamics if the static dielectric permittivity of a 
dielectric plate is finite. 

We complete this section by the consideration of the high-temperature limit. 
Here the zero-frequency term (46) of the Lifshitz formula determines the to- 
tal result. In the configuration of metal-dielectric plates only the transverse 
magnetic mode (for which the metal reflection coefficient is equal to unity) 
contributes to the zero-frequency term. As a result, unlike the case of two 
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metal plates, finite conductivity corrections do not contribute at large separa- 
tions (high temperatures). Thus, for metal and dielectric plates made of real 
materials, equations (47) and (48) obtained for ideal metal and dielectric with 
constant permittivity preserve their validity. 



5. Comparison between analytic and numerical results 



Here we compare the analytic results for the Casimir energy, free energy and 
pressure given by Eqs. (69), (71), (62) and (37) with the results of numerical 
computations using the Lifshitz formulas (14) and the dielectric permittivities 
e^'^{i^i) determined from the tabulated optical data for the complex index of 
refraction. This comparison permits us to find the applicability regions of the 
obtained analytic results for different materials. As an example we consider 
the metal plate made of Au and the dielectric plate made of high-resistivity 
Si. 

The most precise results for (i^i) in the case of Au were obtained in [23] and 
for e^{iii) in the case of Si in [62]. In both cases the data for Im£-^'^(a;) were 
taken from [63] and the dielectric permittivities along the imaginary frequency 
axis were computed by means of the Kramers- Kronig relation, 

oo 



Note that the dielectric permittivity of Si along the imaginary frequency axis 
is equal to its static value {Eq = 11.66) up to the angular frequency of 5 x 
lO^^rad/s and with increase of frequency decreases to unity. The analytical 
results were computed with the plasma frequency of Au equal to cUp — 9.0 eV 
and the characteristic absorption frequency of Si equal to cui — 4.2 eV [63] 
(leV = 1.519 X lO^^rad/s). 

In Fig. 2 we compare the results of analytic and numerical computations of 
the Casimir energy density (A) and pressure (B) at different separations at 
zero temperature. In the vertical axes the quantities 6E = {Ea — En) /En (A) 
and 5Pq = {Po,a — Po,n)/Po,n (B) in percent are plotted where Ea and Po,a 
are the analytic results calculated by Eqs. (69) and (71), respectively, and En, 
Po,n s-re computed numerically using the Lifshitz formula as described above. 
As is seen in Fig. 2, the largest deviations between the analytic and numerical 
results (-4.3% and -7.1% for the energy and pressure, respectively) hold at the 
shortest separation of 100 nm. This is because the plasma model works good 
only at separations larger than the plasma wavelength. At shorter separations 
not some analytic representations for e but the tabulated optical data should 
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be used to obtain precise results. At separations larger than 200 and 300 nm 
\5E\ is less than 0.9% and 0.25%, respectively. As to |5-P|, it is less than 0.9% 
and 0.25% at respective separations larger than 250 and 370 nm. Thus, the 
obtained analytic formulas for the Casimir energy density and pressure at zero 
temperature in between metal and dielectric give rather precise results in a 
wide separation range with a precision at the fraction of a percent. In some 
cases this makes unnecessary much more cumbersome numerical computations 
using the Lifshitz formula and tabulated optical data for the complex index 
of refraction (note that the use of different sets of tabulated optical data also 
leads to about 0.5% differences in the numerically computed Casimir forces 
[17]). 



In Fig. 3 the results of analytical and numerical computations of the relative 
thermal correction to the Casimir energy (A) and pressure (B) at a = 300 nm 
are compared at different temperatures. The relative thermal corrections are 
defined as 



A.F(a, T) ^ T{a, T) - E{a) Ap{a,T) ^ P{a,T) - Poja) 
E{a) E{a) ' Po(a) 



where E{a) and Po{a) are the energy density and pressure calculated nu- 
merically by using the Lifshitz formulas at zero temperature. The analytical 
computations of the thermal corrections are performed using Eqs. (62) and 
(37). Their results are shown by the dashed lines. The numerical computations 
of the thermal corrections are done with the help of the Lifshitz formula at 
zero and nonzero temperatures (solid lines). As is seen in Fig. 3A, the low- 
temperature analytic result for the thermal correction to the energy density 
reproduces the result of numerical computations at T < 20 K. From Fig. 3B 
it follows that the low-temperature analytic expression for the thermal cor- 
rection to the Casimir pressure works good in a wider temperature region 
T < 40 K. The deviations between analytical and numerical results at higher 
temperatures are explained by the fact that in Eqs. (62) and (37) we have 
restrected ourselves by only two and one perturbative orders in small param- 
eter r, respectively. This restriction, however, makes it possible to solve the 
main problem of our interest which has no numerical solution, i.e., to find the 
behavior of the Casimir free energy, entropy and pressure at arbitrarily low 
temperatures. 
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6. Is the Lifshitz formula for configuration of metal and dielectric 
consistent with thermodynamics? 



In Sections 3-5 it was supposed that at zero frequency the dielectric permit- 
tivity of the dielectric plate is finite. It is well known, however, that at nonzero 
temperature dielectrics possess some dc conductivity (Jq = cto(T) which is very 
small in comparison with the conductivity of metals. Usually (see, e.g., [48,64]) 
this conductivity is included into the model of dielectric response by adding a 
Drude-like term in the dielectric permittivity of dielectric, 

i"m=s-(ia + '-^. (75) 



Eq. (75) presents the typical example of the situation discussed in Section 4 
when the dielectric permittivity depends on the temperature as a parameter. 
It can be identically represented in the form 

e^'m = e^'m + (76) 



where P{T) = 2hao{T)/{kBT). The conductivity of dielectrics quickly de- 
creases with temperature, cro(T) ~ cxp{—g/T), where the coefficient g is 
determined by the width of the energy gap A [65]. The magnitude of the 
additional term P{T)/l in Eq. (76) is very smaU. Thus, for SiOs at T = 300 K 
it holds (3 ~ 10~^^ [66]. This makes the role of dielectric dc conductivity 
negligible at alH > 1. 

The question arises on the possible role of dielectric dc conductivity in the 
Casimir interaction between metal and dielectric. As was shown in [45,46,47], 
the Lifshitz formula for the Casimir interaction between two dielectrics cannot 
incorporate the effects of the dc conductivity because then an inconsistency 
with thermodynamics arises. The substitution of Eq. (75) in Eq. (10) leads to 

ff(0,y) = l, ff(0,y)=0, (77) 



instead of Eq. (15). Despite the negligible role of the dc conductivity at all 
I > 1, this could lead to important consequences for the Casimir interaction 
between metal and dielectric. Importantly, the inclusion of the dc conductivity 
of a dielectric plate leads to a discontinuity in the transverse magnetic reffec- 
tion coefficient at zero frequency as is seen from Eqs. (15) and (77). This is 
unlike the case with metals described by the Drude model where the discon- 
tinuity arises in the transverse electric reflection coefficient at zero frequency. 
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To investigate this problem, we substitute the dielectric permittivity (76) in 
Eq. (9) instead of e^{i^i) defined in Eq. (50). For a metal, as in Section 4, 
the dielectric permittivity in Eq. (49) is used. In such a way the Casimir 
free energy jF(a, T) is obtained which takes into account the effects of the 
dielectric dc conductivity. It is convenient to separate the zero-frequency term 
of ^-"(0, T) and subtract and add the zero- frequency term of the free energy 
J^{a,T) calculated with the dielectric permittivity e^{i^i): 



^(a, T) = I ydy [in (l - e"^) - In (l - roe"^) 



00 



+ 





00 



Y: J ydy {in [1 - (0, y)ff{Ci, y)e-y 



+ ln [l-rY{Ci,y)f^iCi,y)e-y]}. 

Here the reflection coefficients f^j_ are found by using Eq. (10) where the 
dielectric permittivities e^{i^i) from Eq. (50) are replaced by e^{i^i) from 
Eq. (76). 

To find the behavior of jF(a, T) at low temperatures, we expand the last in- 
tegral on the right-hand side of Eq. (78) in powers of the small parameter 
j3{T)/l. The zero-order contribution in this expansion together with the sec- 
ond integral on the right-hand side of Eq. (78) are equal to the Casimir free 
energy jF(a,T) calculated with dielectric permittivity e^{iC,i). Calculating ex- 
plicitly the first integral on the right-hand side of Eq. (78), we rearrange this 
equation to the form 

^{a, T) = /-(a, T) - [((3) - Li3(ro)] + Q{a, T) , (79) 



where Q{a,T) contains all powers in the expansion of the last integral on the 
right-hand side of Eq. (78) in the small parameter j3{T)/l equal or higher than 
the first one. The explicit expression for the main, linear in (3{T)/l, term in 
Q{a,T) reads: 



n( ^^ k^T ^ p{T) 7 dyy^e-y 



(80) 
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X < 



(2 - er)a - 2y' 
y2 + Cf{ef> - 1) + ef^y 



ffiCi,y) 
i-ff{Ci:yK{Ci:y)e-y 



^/y' + C?(eP-l)+y 



l-rjf{Ci,y)f1{Ci,y)e-y 



To determine the asymptotic behavior of Eq. (80) when r — > 0, we expand 
the integrated function in powers of r (recall that Q = tI) and consider the 
main contribution in this expansion at r = 0: 



gi(a,T) = 



E 



/?(T) 7 ydye-y 



kBTP{T) 



oo n 



47ra2 ({sEf -l)n 

Performing the summation in / we obtain 

kBTl3{T) 



oo ^—nrl 



E —r- + «^ E e 
1 



7=1 



-nrl 



Qi{a,T) 



Ana' {{e^f - l) 



E 

n=l 



-ln(l-e-"") + 



nr 



(81) 



(82) 



The right-hand side of Eq. (82) can be rearranged with the help of the equahty 



nr 



ln(l - e-"0 + — — - = -lnT + l-lnn + O(t^). 



(83) 



As a result it holds 

Qi(a,T) = 'it't'^ . T(3(T)\nr + T(3(T)0(r'). 



2na' {{e^f - l) 



(84) 



Taking into account that I3{T) ~ (l/T") exp(— g'/T), we arrive at 

Qi{a,T) r^e-^'^hiT. (85) 



From Eq. (85) it follows that both Qi{a,T) and its derivative with respect 
to T go to zero when T vanishes. The terms of higher orders in the small 
parameters P{T)/l and r omitted in our analysis, go to zero even faster than 
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Qi. Thus, the quantity Q{a,T) in Eq. (79) and its derivative with respect to 
T have zero hmits when the temperature goes to zero. 



Now we are in a position to find the asymptotic behavior of the entropy in the 
configuration metal-dielectric with included dc conductivity of the dielectric 
plate. Using Eq. (12), we obtain from Eq. (79) 

S{a, T) = S{a, T) + ^ [C(3) - Li3(ro)] - (86) 



where S{a,T) is defined in Eq. (72). In the hmit T ^ Eq. (86) results in 
'5(a,r) = ^[C(3)-Li3(ro)]>0. (87) 



This equation implies that in the configuration of metal-dielectric with in- 
cluded dc conductivity of the dielectric plate the Nernst heat theorem is vi- 
olated. Previously the analogous result was obtained [45,46,47] for the con- 
figuration of two dielectric plates with frequency- dependent dielectric permit- 
tivities. It is easily seen, that Eq. (87) is preserved, if, instead of the plasma 
dielectric model (49) , the Drude dielectric function is used in the case of metal 
plate with a perfect crystal lattice. If the metal plate has impurities, the an- 
alytical derivation cannot be performed, but numerical computations lead to 
the same positive value of the entropy as in Eq. (87). Thus, the Lifshitz theory 
becomes inconsistent with thermodynamics when the dc conductivity of a di- 
electric plate is taken into account. This suggests that the actual low- frequency 
behavior of the dielectric properties is not related to the phenomena of van 
der Waals and Casimir forces and should not be included into the model of 
dielectric response. 

Recently [28] this theoretical conclusion was confirmed experimentally in the 
measurement of the difference Casimir force acting between Au-coated sphere 
and Si plate illuminated by laser pulses. The difference of the Casimir forces in 
the presence and in the absence of pulse was measured using an atomic force 
microscope. In the absence of laser pulse the concentration of charge carriers 
was of about 5 x lO-'^^cm"^ (higher-resistivity Si), but in the presence of pulse 
this concentration has been enhanced up to 2 x 10^^ cm^'^ (lower-resistivity 
Si). The experimental data were compared with two theories. The first theory 
used an assumption that in the absence of laser light Si possesses a finite static 
dielectric permittivity — 11.66 (see Section 5). The second theory took into 
account the dc conductivity of Si in the absence of laser light like it was done 
above in Section 6. The first theory was found to be in excellent agreement 
with data, whereas the second theory was excluded at the 95% confidence 
level within the separation region from 100 to 200 nm [28] . Thus, the inclusion 
of the dc conductivity of dielectrics and high-resistivity semiconductors in the 



28 



model of dielectric response is not only inconsistent thermo dynamically but 
is also in contradiction with experiment. 



7. Conclusions and discussion 

In this paper we have obtained the analytic expressions for the Casimir free 
energy, pressure and entropy at low temperatures in the configuration of one 
metal and one dielectric plate. Different models of the dielectric response for 
both metal and dielectric were considered: the simplified model of an ideal 
metal and dielectric with constant dielectric permittivity, and the realistic 
model of a metal described by the plasma model and a dielectric described in 
the Ninham-Parsegian representation for e. To derive the asymptotic expres- 
sions at low temperatures, the perturbation theory in the small parameter r 
was developed which is proportional to the product of separation distance and 
the temperature. The analytic expressions for the main physical quantities in 
the limit of high temperatures and at zero temperature were also obtained. 
The analytic results were compared with numerical computations using the 
Lifshitz formula and tabulated optical data, and good agreement was found. 

The fundamental conclusion arrived in the paper is that the Lifshitz theory 
applied to the configuration of one metal and one dielectric plate is in agree- 
ment with thermodynamics if the dielectric permittivity of a dielectric plate at 
zero frequency is finite. In particular, it was shown that the Casimir entropy 
goes to zero when the temperature vanishes, i.e., the Nernst heat theorem is 
satisfied. This conclusion cannot be reached numerically, because it is impos- 
sible to perform numerical computations at arbitrarily low temperature and 
their precision is always restricted. On the contrary, it was shown that, if the 
dielectric permittivity of the dielectric plate at zero frequency turns to infinity 
(i.e., small dc conductivity of a dielectric material is taken into account), this 
leads to a nonzero value of the Casimir entropy at zero temperature, i.e., to 
a violation of the Nernst heat theorem. The inclusion of the dc conductivity 
of high-resistivity semiconductors in the model of dielectric response was also 
recently shown to be in contradiction with experiment [28] . What this means 
is that to avoid contradictions with thermodynamics and experiment, one 
should not include the actual conductivity properties of dielectric materials 
at very low, quasistatic, frequencies into the model of the dielectric response 
(this phenomenological prescription was obtained previously in [45,46,47] for 
the case of two dielectric plates made of real materials) . 

The above conclusions can be discussed in the context of the formalism of ther- 
mal quantum field theory in Matsubara formulation where the zero-frequency 
term plays a separated role and calls for an adequate interpretation. It is 
common knowledge that the zero-point energy of quantized fields contains os- 
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cillations of any frequency. It would be hard, however, to imagine the presence 
of a zero- frequency (i.e., constant) field in the vacuum state. This suggests that 
the zero-frequency term in the Matsubara summation could be understood not 
literally but as a mathematical limit to zero from the region of much higher 
frequencies [the characteristic frequency of the Casimir effect c/(2a), and the 
thermal frequency kBT/h] which determine the physical phenomena of the van 
der Waals and Casimir forces. As was discussed below Eq. (76), in the region of 
characteristic and thermal frequencies the effects of dc conductivity contribute 
twelve orders of magnitude less than e^. Because of this, the dc conductivity 
may be considered as not related to the Casimir forces and be not included 
into the zero-frequency term of the Lifshitz formula. This phenomenological 
prescription is not the fundamental resolution of the problem, which remains 
unknown, but following it one avoids contradictions to thermodynamics and 
experiment. 

Note that the problems discussed above for the configuration of metal-dielectric 
are of different nature than those arising for two metals. For metals the con- 
centration of charge carriers only slightly depends on the temperature. The 
validity of the Nernst heat theorem in the Casimir interaction between two 
metals is caused by the scattering processes of free charge carriers on phonons, 
impurities etc. For the perfect Drude metals with no impurities, relaxation 
goes to zero when the temperature vanishes and the Nernst heat theorem is 
violated [33,43,44]. On the contrary, for dielectrics the concentration of charge 
carriers quickly decreases to zero when the temperature vanishes. Here the 
violation of the Nernst heat theorem does not depend on the scattering pro- 
cesses and is caused by the inclusion of the infinitely large dc conductivity. 
In the formalism, for two metals a discontinuity in the refiection coefficient of 
the transverse electric mode at zero frequency arises, whereas for a metal and 
dielectric a discontinuity holds in the reflection coefficient of the transverse 
magnetic mode. These differences are reflected in the fact that even the sign 
of the entropy at T = for two perfect Drude metals and for one metal and 
one dielectric with included dc conductivity are opposite (negative for two 
metals and positive for a metal and dielectric). 

The obtained results are topical for the interpretation of recent measure- 
ments of the Casimir force between metal sphere and semiconductor plate 
[25,26,27,28]. Semiconductors suggest a wide variety of the conductivity prop- 
erties ranging from metallic to dielectric ones. In the application of the Lif- 
shitz theory to the metal-semiconductor test bodies the model of the dielectric 
response for a semiconductor should be chosen to satisfy the Nernst heat the- 
orem and other fundamental physical principles. This can be done by using 
the proposed phenomenological prescription. A more fundamental resolution 
of the discussed problems may go beyond the scope of the Lifshitz theory. 
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Fig. 1. The correction factor ipDM to the Casimir energy density (solid line) as a 
function of the static dielectric permittivity. The long-dashed lines 1, 2 and the 
short-dashed line show the coefficients Ci, C2 and B, respectively, in Eq. (69) for 
the Casimir energy density between plates made of real metal and dielectric. 



35 




Fig. 2. The relative differences between the analytic and numerical results for the 
Casimir energy density (A) and pressure (B) at zero temperature versus separation. 
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Fig. 3. The relative thermal correction to the Casimir energy density (A) and pres- 
sure (B) versus temperature at a separation a = 300 nm. Solid lines show the results 
of numerical computations and the dashed lines are obtained using the analytic 
asymptotic expressions. 
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